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O . ABSTRACT 

' The effects of radiation drag force on the structure of relativistic electron-positron 
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and electron-proton outflows are considered within the one-fluid approximation for 



' quasi-monopole cold outflow. It is shown that for a Poynting-dominated case the drag 

, force does not change the particle energy inside the fast magnetosonic surface. In this 

'"^ ' region the action of the drag results in a diminishing of the Poynting flux, not the 

particle flux. Outside the fast magnetosonic surface, for intermediate photon density 
the drag force may result in additional acceleration of the plasma. This acceleration 
is a result of the disturbance of magnetic surfaces under the action of the drag. At 
' even larger distances particles are not frozen into the magnetic field and the drag force 

, decelerates them efficiently. 

' In the case of extreme photon densities, the disturbance of magnetic surfaces be- 

. comes large and the drag force changes the total energy flux significantly, the particles 

' becoming nonrelativistic. 

. We find that for Active Galactic Nuclei the photon density is too low to disturb 

the parameters of an ideal MHD outflow. The drag action may result in additional 
^ , acceleration of outgoing plasma only for central engines with very high luminosities. 

For cosmological gamma-ray bursts the drag force can strongly affect the process of 
formation of a Poynting-dominated outflow. 



^ ' Key words: acceleration of particles - MHD - galaxies: active - gamma-rays: bursts 



1 INTRODUCTION 

. . 1 Magnetohydrodynamic (MHD) models are now developed intensively in theories of the magnetospheres of rotating supermas- 
sive black holes (M ~ 10^-10^ Mq, Bq ~ 10* G), which are believed to reside in central engines of Active Galactic Nuclei 
(AGNs) and quasars (Begelman, Blandford & Rees 1984; Thorne, Price & Macdonald 1986). In particular, it is the MHD 
model that is the most promising in the problem of the origin and stability of jets. Indeed, the MHD approach explains both 
the energy source (the rotational energy of the compact object) and the mechanism of the energy and angular momentum loss 
(for an overview, see e.g. Blandford 2002). Observational evidence in favor of MHD models was recently found in the possible 
presence of toroidal magnetic fields in jets (Gabuzda et al 1992; Gabuzda et al 1999). Magnetically-dominated outflows are 
also believed to be responsible for the energy transport in cosmological gamma-ray bursts (Meszaros & Rees 1997; Lee, Wijers 
& Brown 1998; van Putten & Levinson 2003), when energy is released in the merging of black holes or neutron stars (M ~ Mq, 
Bo ~ lO^'^ G). 

It has been suggested that the density of photons in the vicinity of the central engine is so high that they may drastically 
change the characteristics of the ideal MHD outflow. For example, they may result in extensive e^e~ pair creation (Svensson 
1984), acceleration of low-energy pairs by the radiation drag force (Phinney 1982; TuroUa, Nobili & Calvani 1986; Beloborodov 
1999) and deceleration of high-energy particles (Melia & Konigl 1989; Sikora et al 1996). In other words, a self-consistent 
consideration should take the drag force into account. 

So far the two processes - the ideal MHD acceleration and the action of external photons - have been considered 
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independently. The first step to combine them was made by Li, Begelman & Chiueh (1992). In particular, it was shown 
how the equations can be integrated in a conical geometry (which is impossible in the general case). On the other hand, 
the consideration was performed within the approximation of a fixed poloidal magnetic field. Under this assumption the fast 
magnetosonic surface of a cold flow is shifted to infinity (Michel 1969; Kennel, Fujimura & Okamoto 1976; Lery et al 1998). 
As a result, it was impossible to analyze the effects of radiation drag on the position of the fast magnetosonic surface and the 
properties of the outflow outside this surface. 

The main goal of this paper is to determine more carefully the radiation drag effects on a magnetically-dominated outflow. 
To describe analytically the effects of radiation drag, including simultaneously the disturbance of the magnetic surfaces we 
consider here a quasi-monopole outflow. For AGNs, such geometry in the immediate vicinity of the central engine was recently 
confirmed by direct observations (Junior, Biretta & Livio 1999). In other words, in the zeroth approximation (i.e., without 
drag) we use the analytical solution for a magnetically-dominated MHD outflow (Beskin, Kuznetsova & Rafikov 1998, hereafter 
Paper I), in which the fast magnetosonic surface is located at a finite distance from the origin. 

For simplicity we consider the following model of the radiation field in the vicinity of the central engine. First, we notice 
that for ultra-relativistic particles the energy of a photon propagating nearly along the particle trajectory remains almost the 
same after a collision. This means that the drag force from these photons is small. Thus, only the isotropic component of 
the photon field contributes substantially to the drag force. Hence, in our geometry with a strong central source of photons 
and a monopole outflow of particles, only a small fraction of photons (the isotropic component of the photon field) interacts 
efficiently with the particles, producing inverse Compton photons with energies 8ic ~ 7^fph. 

The isotropic component can be produced, firstly, by the outer part of the accretion disk and, secondly, by external 
sources. It can be modeled as 

U =Uiso = UA[^ " + C/ext, (1) 

where J?l = c/O is the radius of the light cylinder, Ua = U{R\.) = Ltot /(47r i?£c), and n w 3 (for more details see, e.g., 
Sikora ct al 1996). Here the first term describes the radiation from the outer parts of the disk, Trad > r, while the second one 
corresponds to the homogeneous external radiation. For AGNs this can be due to clouds located at a distance rdoud ~ Ipc 
from the central engine and reradiating kL^ot of the total luminosity {k ~ 10%). In this case 

U...=k^^. (2) 

^" ' cloud'- 

However, this model only makes physical sense at distances less than rdoud, and the term vanishes at larger distances. 

Finally, as some arguments exist both in favor of (Reynolds et al 1996; Hirotani et al 1999) and against (Sikora & 
Madcjski 2000) the leading role of e+e" plasma in relativistic jets, in what follows we consider both electron-positron and 
electron-proton outfiows. 

In Section 2 we formulate the basic equations describing a quasi-monopole outflow of relativistic plasma in two-fluid and 

one-fluid approximations. Then in Section 3 we analyze the main properties of an electron-positron outflow. A similar analysis 
for electron-proton plasma is produced in Section 4. Finally in Section 5 we consider the effects of radiation drag for real 
astrophysical objects. 



2 BASIC EQUATIONS 

2.1 The Two-Fluid Description 

We consider an axisymmetric, stationary outflow of two-component cold plasma from the magnetosphere of a rotating body 
with a split monopole poloidal magnetic field. This geometry can be realized in the presence of an accretion disk separating 
the ingoing and outgoing magnetic fluxes (Blandford & Znajek 1977). In the hydrodynamic approximation, the structure 
of the flow is described by Maxwell's equations and the separate equations of motion for positively and negatively charged 
particles: 

VE = 4npe, V xE = 0, 

VB = 0, V X B = —j, (3) 

c 

(t;± V)p± = ±e(^E + ^xB^+ F^^^. 

Here E and B are the electric and magnetic fields, pe ~ e{n'^ — n^) and j = e{n'^v'^ — n~v~) are the charge and current 
densities, and and are the velocities and momenta of the charged particles. Finally, F^^^^ are the radiation drag forces 
which, for an isotropic photon field, have the form 
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4 V 



± 



m± 



(4) 



-^drag 3|^±| 

Here 7^ are the Lorontz factors, rri-^- is the mass of the positively charged particles (positrons or protons) and m_ = rric is 
the mass of the negatively charged particles, ut = (87r/3)(e^/meC^)^ is the Thompson cross section. To close system (3) the 
continuity equations V{n^v^) = should be added. It is enough to add the equation for one component, e.g. 

V{n+v+) = 0, (5) 

because the continuity equation for the second component then follows from (5) and the current continuity equation Vj = 
(which, in turn, results from Maxwell's equation V x B = {An/c)j). 
In the limit of infinite particle energy. 



(0) 



7 = oo, t;;"' = c, = 0, V, 

so that 

Pe= Ps^ COS 0, jr = PsC—^ COS 6, je = = 0, 

the monopole poloidal magnetic field 



(0) 



0, 



B 



(0) 



is an exact solution of Maxwell's equations. In this case. 



= E. 



(0) 



-Bo smf 

c r 



0, 



(6) 
(7) 

(8) 

(9) 
(10) 



which is just the well-known Michel (1973) solution. Hero Bo and ps are the magnetic field and charge density on the surface 
r = R <^ Rh, and the angular velocity is f2 = 2nc\ps\/Bo. The limit 7 — > oo corresponds to zero particle mass in the force-free 
approximation. 

It is also convenient to introduce the electric field potential ^{r,9), so that E = — V$ and 
Q.R^Bo 



(0) 



■ cos 6, 



(11) 



and the flux function ^{r,9), so that the poloidal magnetic field 
„ V* X 



2nr sin 9 

and = 27rB(,i?2(l - cosO). 

Then the dimensionless corrections T;*(r, 
the following form: 

+ fiBo R^ 



(12) 



^*(r, 0), 5{r,6), £f{r,9) and Ci'i'y^) for the case v 7^ c can be introduced in 



^{r,9) 
^{r,9) 

Br 

Be 

Bip 

Er 
Eg 



2nce 
OBo R^ 



A - ^ cos 9 + r]+{r, 9) 



A -I- - cos 6* -f Tj 



27rce 

c[l-e±(r,e)] , vf =c(,^{r,9), 
^^^-^[-cose + 5(r,6>)], 



2'kBoR^ [1 - cos 61 + e/(r, 9)] , and thus 



s^Rl I e df 

sin 6* 89 
BoR^ df 

— £ — . 

r sin Or ' 
Bo — -[- sm9-ar,9)]: 



nBpR^ 85 
c 9r' 
ClR^Bo 



(13) 

(14) 
(15) 

(16) 
(17) 
(18) 

(19) 

(20) 

(21) 

(22) 
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Here A » 1 is the multiplication parameter (A = ens/\pa\, where ng is the number density of particles on the surface r = R). 
For A < 1 the approach under consideration is not valid. In what follows wo consider for simplicity the case A = const. Such 
a choice corresponds to a constant particle-to-magnetic flux ratio k = const. 

Switching to dimensionless variables, we use below the dimensionless radius x = r/Ri, = ril/c and the dimensionless 
drag force: 

4 axC/iso , 



3 ClnieC ^ ' ' ' ^^^^ 
Now, substituting (13)-(22) into (3) we obtain, to first order in all the correcting functions, the following system of equations: 



d 



iCsme) = 2{r,+ -r,-) 



(a- i cose) (^A+icose) e 

.dS 



V ) + ■ 



dx 



1 d 
sin e de 



,dS 



= 



sm a 

X 



dx'^ 



±1. 

^2 an 



sin 6 





2 


dx 


X 


df\ 


2 ■ 


86 ) 


X . 



A — - cos 



A -|- - cos ( 



Tx (^«"^") + 



dx 



rrip J \ rup J \ x off x x x^ ) 

sin6» \ 



d_ 

dx 



d-x («-"^") + 



rup 



^1 +4Aa 



- +4Aa — 
Up J V^p/ ^ 

(7-) =-Fd--4Aa( 

rric \ / 
rUp J V 

-4Afr 



(- 



d5_ _ 
dx 
d5_ _ 
dx 

1 df 
X sin dx 

1 df 
X sin 6 dx 



Here 



QeBoR^ 



> 1 



(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 

(34) 



is the Michel (1969) magnetization parameter describing the particle-to-electromagnetic energy flux ratio W^part/M^em = 
{mp/me)'y/cr. Hence, for a Poynting-dominated flow we have 7 ^ cr. As we see, the disturbances of the particle density r?"*" 
and ri~ enter equations (24)-(33) only in the combination rj'^ —ri~ . Therefore, the system (24) -(33) is closed, but the equation 
of mass continuity (5) is necessary to determine r]'^ and ri~ separately. The system (24)-(33) differs from the one considered 
by Beskin & Rafikov (2000, hereafter Paper H) only by the additional drag terms in the r.h.s. of (28)-(33). 



2.2 The One-Fluid Limit 

Equations (24)-(33) describe the flow in the two-fluid approximation. We now reduce the complete system of equations (24)- 
(33) to consider the one-fluid approximation. In the case of an electron-proton outflow there is a small parameter nie/mp ~ 
10~^ which allows us to neglect the electron mass and thus to proceed in a standard way to the one-fluid approximation. 
But as demonstrated in Paper II, this can be done in the electron-positron case for a magnetically-dominated (cr ^ 1) dense 
(A ^ 1) plasma as well. Because the nonhydrodynamic components of the velocity are small in this case (cf. Melatos & 
Melrose 1996), 

^~A-V-^/^ ^^X-\ ^^x-\-^'\ (35) 

we can set = (~ = {i = r,9, ip), where is the hydrodynamic velocity. In this limit we also have 

- A-V-^/^ (36) 

As a result, for a magnetically-dominated outflow with cr ^ 1 and A ^ 1 in the one-fiuid approximation, 

S = ef. (37) 
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Finally, equations (28) and (29) together with (37) in the limit Au ^ 1 under consideration give another useful one-fluid 

relation: 

X 39 X X 

this equation being the same as in the drag-free case. On the other hand, as demonstrated in Paper II, 

|^W(T-'/^ (39) 

Hence, in the one-fluid approximation one can set = so that 



7 = T2- (40) 



.2 ^ 1_ 

2^T- ~ 

It is necessary to stress that in some cases the monopolc geometry allows one to consider separately the set of equations 
describing the particle energy and the set of equations resulting in the Grad-Shafranov (GS) equation, which determines the 
disturbance of magnetic surfaces. Thus we can consider particle energy without formulating here the general form of the GS 
equation. Some asymptotics of the GS equation are discussed below. 



3 THE ELECTRON-POSITRON OUTFLOW 

3.1 Integrals of motion 

In this section we consider the properties of the electron-positron outflow when rrip = nie = m. In our simple geometry with a 
(split) monopole poloidal magnetic field in the zeroth approximation the particle motion can be considered as radial. It allows 
us to integrate some equations along the r coordinate axis. This feature was first demonstrated by Li et al (1992). 
Indeed, combining (30) and (31) with (26) one can obtain in the one-fluid approximation 

^ = Te^^e^-^e^^-^-^-^e'^ /^(.V(.')d.'. (4i) 

Here 7in is the Lorentz factor near the origin, xq = QR/c, and the integration constant l{9) describes the disturbance of 

the electric current I(R,0) = Ia [sin^ 9 + /(6)] on the surface r = R. The integration constant must be determined from the 
critical conditions on singular surfaces. The compactness parameter Ia is deflned by the relation 

/a = 1^1^ (42) 
3 mecQ ' 

where Ua = U{Ry,) is again the first term in (1) on the light cylinder r = i?L. Finally, we also denote 

^""-'-s^^' ^^^^ 

so that (1) for a; > 1 can be rewritten as 

u{x)=x-" + le.t/lA. (44) 

For a; < 1 it is natural to set u{x) = 1. 

Expression (41) is just another way of writing the diminishing of the Bernoulli integral Eb along a magnetic field line 
as a result of the drag force. Within the Grad-Shafranov approach the full energy loss is determined by = J Esd'^f. Thus 
relation (41) can be rewritten as 

An f 

Esir) = Eb{R) - — / Fdragdr, (45) 



where the energy flux per unit magnetic flux Eb has the standard form Eb ~ 0.1/2110 + jmc^K. Here / is the total electric 
current inside the magnetic tube and n is the particle-to-magnetic flux ratio (see e.g. Beskin 1997 for details). 
In what follows we consider the case 7in ~ 1, i.e. 

Tin < ff, (46) 

when additional acceleration of particles inside the fast magnetosonic surface takes place (see e.g. Paper I). In the case 7;^ » a 
corresponding to ordinary pulsars, the particle energy remains constant (7 = 7in) on any way up to the fast magnetosonic 
surface (Bogovalov 1997). 

The other two integrals of motion, namely the conservation of angular momentum separately for electrons and positrons. 
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can be obtained from equations (30)-(33): 

7+(l-a;sine^J) =7+ -/a / u{x'){l - x' sinei+){'i+fdx' - A\a{5 - ef), 



7 (1 - xsin^^^ ) = 7in - Za / M(a;')(l - a;' sin^^^ )(7 )'^Ax' -\- AXa{5 - ef). 
In the one-fluid approximation this gives 

7(1 — xsin^^^) = 7in — /a / u{x')(l — x' sm9^,p)^^{x')dx' . 

J Xg 

At large distances, where 7 ^ 7in, we have 



(47) 
(48) 

(49) 
(50) 



It can be seen from (49) that the presence of the drag force (the term proportional to ^a) makes this approximation more 
accurate. 



3.2 Fast Magnetosonic Surface 

Substituting ^ from (41), from (40), and from (49) into (38) we get the following equation to determine the position of 
the fast magnetosonic surface r = rp: 

^2" 



ad tand smd asiad 



1 
2? 



+ 



(7 - 7in + Jl) 

2x272 sin^ e 



7 - 7in + Jl 

a;27 sin ^ 



0. 



Here 
J = 



/a 



'u(x-')7"(x-')da;'. 



rx 

Jl = /a / w(a;')7^(a;')(l — a;'sin0£<^)da;'. 

J xn 



(51) 



(52) 



(53) 



But, in accordance with (49) and thereafter, the terms Ji and 7in can be neglected far from the origin of the flow r S> i?L- 
Also, since l{0) ^ a~*^^ (see Paper I), the term with l{6) can be omitted as well. Hence, one can write down the following 

algebraic equation for the Lorentz factor 7: 



7 

where 

P = - 



df 



J 



2ef cos 6 — £ sin 



(54) 



(55) 



This differs from the drag-free case by the additional term J/a. 

Equation (54) allows us to determine the position of the fast magnetosonic surface r = ry and the energy of particles on 
this surface 7f = 7(rp). Indeed, as the fast magnetosonic surface is the X-point, we flnd the exact solution for coinciding 
roots 



1/3 • 2/3 , 
7p = (7 ' sm ' ( 



(56) 



which does not depend at all on P, or thus on the drag force. On the other hand, both the numerator and the denominator 
of the derivative d7/da; must be equal to zero on the surface r = rp: 

d7 _ 7^ (a:dP/da: - x~^) 

^d^ ~ 37 - CT (2P a;-2) ■ ^ ^ 

It gives 

{P + x-^)f ~ a-''^'\ |P|f « Ix-'If. (58) 
Now using (55), one can see that the conditions for weak and strong drag are, respectively, I a <IC ^cr and I a ^ lev, where 

/cr = a^/^ (59) 

In other words, for Ia <Si Icr the flow remains the same as in the drag-free case, and 

XF « (7'''^sin"^/^6i, (60) 
(£/)p ^ o-^l\ (61) 
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On the other hand, for a high enough photon density Ia S> la one can obtain, for n = 3, 

(sfh « (63) 

Thus, as we can see, the energy jFmc^ do not depend on the drag. The disturbance of the magnetic surfaces increases with 
increasing I a, but remains small for Za "C cr. 

For a very dense photon field, according to (41) and (63), the disturbance of the magnetic surfaces becomes of order unity 
for I A ~ Zmax where 

^max = a, (64) 

i.e., when the disturbance of the Bernoulli integral AEb = Esir) — Eb{Rl,) (45) is on the order of the total energy flux Eb, 
or in other words when the drag force substantially diminishes the total energy flux of the flow. But this means that the 
disturbance of the magnetic surfaces becomes large when the energies of particles in the vicinity of the fast magnetosonic 
surface become nonrelativistic. This feature is well known for a drag-free flow both within numerical (Sakurai 1985; Bogovalov 
1997) and analytical (Bogovalov 1992; Tomimatsu 1994; Paper I) considerations. Of course, our analysis within the small 
disturbance approach can only demonstrate the tendency. For this reason, within our approach a = const. In reality, since 
a = EB/{mc^K), a decrease of the energy flux results in a decrease of a as well. 

Finally, it is necessary to stress that the results depend signiflcantly on the assumed model of the isotropic photon density. 
In particular, for the general power-law dependence of the photon density for n < 3 in (44), one can obtain for the position 
of the fast magnetosonic point and the disturbance of the magnetic surfaces 
/ ; \ 1/(5-™) 

« cr^/^(^) , (65) 

/;, \2/(5-n) 

(£/)f « (^) (66) 

instead of (62) and (63). This takes place for I a > la, where la is the photon density capable of disturbing the magnetically- 
dominated outflow: 

Z„ = (67) 

A large disturbance (e/)F ~ 1 can only be realized for a very high photon density I a > /max, where again 

/max = a. (68) 

For n > 3 (when the value J defined by (52) is determined by the lower limit of integration and does not depend on n) we 
have (e/)F Ia/ct, xf « (ct/Za)^^^, Ict = c^^^, and Zmax = cr, i.e. same as for n = 3. According to (41), (e/)F ~ 1 occurs 
when AEb/Eb ~ 1 for all values of n. 



3.3 The Structure of the Flow at Small and Lcirge Distances 

3.3.1 Inner region 

We first consider the structure of the flow well within the fast magnetosonic surface r <C rr- As one can see from (54), for 
r < rp 

^ « (69) 

In particular, for Ia ^ Ict one can neglect the term P in the denominator. On the other hand, for Ia ^ la in the immediate 
vicinity of the fast magnetosonic surface the negative term —J/a in P should to be included into consideration. As a result, 
the particle energy increases abruptly up to the value ~ cr^^^. Nevertheless, for r <C rp the value P can be omitted for 
arbitrary compactness parameter I a- It also means that for r rp it is possible to neglect the flrst two terms in (38). Thus, 
in the internal region r <C rp we have 

& = (70) 
a; sm ^ 

Now using relations (50) and (70), one can obtain 

7^ = -|- x^ sin^ 9 x^ sin^ 0, (71) 



c - \/7il + sin^ 8 - _ 1 

xsin^^l^ + x^ sin^e xsm.e 
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^ Von. + -';^'"'^- Tin ^ 1 

a;2 sin^ 6l^7?„ + sin^ 6> a;2 sin^ 61 ' 

in full agroomont with the ideal MHD approximation. Hence, one can conclude that in the internal region r < rp the radiation 
drag does not affect the particle motion. Here the universal value of the Lorentz factor (71) corresponds to the drift velocity, 
and it depends on no external disturbances. Indeed, using the frozen-in condition E + vxB/c = 0, we obtain for the drift 
velocity 

In our case, however, according to (16)-(22), we have 

« E^, El « E'^/{x^ sin 9). (75) 

These relations immediately lead to the previously-stated asymptotic behavior (71). In particular, this means that disturbance 
of the magnetic surfaces plays no role in the determination of the particle energy. For this reason, it is possible not to consider 
the radiation drag corrections to the field structure for r <C rp- 



3.3.2 Outer region 

In the other limit, well outside the fast magnetosonic surface (r » rp) equations (41), (38), (24)-(25), and (40) can be 
rewritten in the form 



J xn 



(76) 



C = e^ + sme^r, (77) 



Here wo neglect the small values l{0), 7in/7 and ^^p, and set 5 = ef. As we can see, the only correction compared to the 
drag- free case is the last additional term in (76). Without the drag force the system (76)-(79) results in a very slow increase 
of the particle energy (7 oc In^/^ r) and actually in the absence of collimation, sf ~ cr~^^^ In^^^ r (Tomimatsu 1994; Paper I). 
Unfortunately, because of the nonlinoarity of the system (76)-(79), in the general case it is impossible to reduce it to the GS 
equation for the magnetic disturbance ef only. For this reason, in what follows we present only an asymptotic representation 
of the GS equation. 

a) Low photon density I a ^ la- 

In this case the action of drag results in a small correction to the particle energy. Only at very large distances, where 
the first term in (44) can be neglected, does the drag significantly change the energy of particles. Indeed, for the decreasing 
component of the isotropic photon density {u{x) ~ x~" in (44)) and for an almost constant particle energy the drag term 
in (76) does not increase with the distance r. Here the increase of the drag term is due to the homogeneous part of the 
photon density (Zext/^A in (44)). Clearly, this can be realized at large enough distances (r > Ri,{lA/le^t)o'^^~"^^^), where the 
contribution from external photons is the leading one. Thus, we see that the action of the drag force can be significant for a 
high enough density of external isotropic photons in the vicinity of the compact object. 

Now we consider this asymptotic region in more detail. Neglecting the term containing 7 in (76) and the first term in 
(44), one can obtain for the generalized GS equation 



d / 2d.f\ 1 d 
Stt- \ X 



o ■ 2n Of ■ 3n9f Sm"^ 6 f 2/ /•,, / 

2esm tfcosoj—ssm 0— text / 7 {x )ax 



= 0. (80) 



dx \ dx , 

One can seek the solution of this equation in the form 

ef[x,e) OCX"' , -i{x,e) <xx°'-' . (81) 
Substituting these expressions into (80) we obtain 

ef{x,e) = fcl(e)4^^-l/^xl/^ (82) 
7(x,e) = fc2(e)C/V^/^x-^/^ (83) 

where k\(9) ^ k2{9) ~ 1 describe the 9 dependence. This takes place for r > rh, where 
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Figure 1. The radial dependence of tlie Lorentz factor ^{r) for low (a) and high (6) photon densities for n = 3. The dashed lines 
correspond to the drag-free case. 



^ i?LCxt^"'/'- (84) 



For rp < r < rh the action of drag is negligible. 

As we see, for r > rh the drag force results in a decrease of the particle energy and an additional coUimation of the magnetic 
surfaces outside the fast magnetosonic surface. These asymptotic solutions, however, do not seem to be realized in reality. 
Indeed, as one can check, rh ~ (10® — 10^)iiL > r-cioud- But for r > rcioud the model that we adopted for the homogeneous 
component of the photon density is no longer valid and thus the asymptotic solutions (82) -(83) do not apply. This conclusion 
is made under the assumption that the homogeneous component is produced by the reradiating clouds surrounding the central 
engine. One could imagine producing a homogeneous component of the photon density in the host galaxy. However, in this 
case the photon density ?7ext is too low and thus the collimation distance rh is too large for our approach to be valid in the 
region r rh. 

The characteristic radial dependence of the particle energy in the presence of the drag is demonstrated in Fig. la. One 

can conclude that for a low photon density the action of the drag force is very weak unless photon field is present out to 
distances r ~ rh. In the latter case the drag force efficiently reduces the particle energy beyond r ~ rh. 

It is necessary to stress here an important property of the fast magnetosonic surface for I a < Ict- Introducing a small 
disturbance of the particle energy resulting from the drag 

7 = 70 + A7, (85) 
where 70 is the Lorentz factor of particles without the drag force, one can obtain from (76)-(79) 

-sin^6l^+e^-2ecos6iA/+iA7=-. (86) 
% 06 a a 

As equation (78) results in eAf ~ A7/70, we see that outside the fast magnetosonic surface (where 70 > a^^^) 
mc'A'yKi- j Fdragdr. (87) 
On the other hand, within the fast magnetosonic surface we have 

mc^ I A7I < y Fd.agdr, (88) 

in full agreement with (71). This means that inside the fast magnetosonic surface the decrease of the total energy flux Eb 

results from the decrease of the electromagnetic energy flux rather than the particle flux. Outside the fast surface the photons 
act upon the particle directly. In other words, for Za < lev the fast magnetosonic surface separates the regions in space where 
particles are strongly or weakly frozen into the electromagnetic field, 
b) High photon density I a > la- 

For a high photon density Ia > Ict the drag force significantly changes the energy of the outgoing particles. As is shown 
in Fig. 16, the increase of the particle energy continues outside the fast magnetosonic surface up to the maximum value 

- ^ 4(n-2)/(5n-n^) 



In 



/ 7 \ '±{11. — — 11 I 

For n < 2 there is no additional acceleration. Only at larger distances (r > rcr) does the radial dependence of the energy of 
the outgoing plasma become similar to the previous case. 
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Figure 2. Saturation energy 7max»7ic^ determined by a numerical integration of (93) for n = 3 and n = 4 (points). The lines correspond 
to the analytical estimate (89) For i a = 'max , n > 3 we have 7sup = c^''^ • 



Indeed, outside the fast magnetosonic surface one can neglect the term sin^^^ in (77). As a result, the system (76)-(79) 
can be rewritten in the form 

'%-'7^f = ^^-^^^a/ u{x')^\x')dx', (90) 

89 tan 9 asm9 asm9 I 

J XQ 



=0 (91) 
dx V dx) 2smHde [ ' ^ ' 



which results in 

dx \ dxj dx ^ ' 



„ „9 /sin^eX . .8^ 
zcosp— — ^ — smf 



72 / 89"^ \ 72 



(92) 



For Ia > Icr the Lorentz factor 7 increases with x, and the r.h.s. in (92) can be omitted. This gives 

:r2|^+«Ax2-V = C(0). (93) 

Here the integration constant C{6) ^ Iax^^^y for 'a S> lev actually does not depend on the boundary conditions d7/da; at 
X = xf- In what follows we do not consider the ^-dependence. 

Now analyzing equation (93), one can find that for a; <C Scr, where 



/7a \ 8/(Bn-n2) 

Xcr , (94) 

the particle energy increases as 

/ „. \ (n-2)/2 

7^^'^'''' — (95) 

\Xy J 

for arbitrary n. As to the disturbance of magnetic surfaces e/, it remains approximately constant: ef « (/a/o")^^'^""'. 
On the other hand, in the saturation region x ^ Xct, where one can now neglect the term iAa;^~"7^ in (93), we have 

7(x) = C (-1 - i) , (96) 

\X\, xj 

where .Tb ~ Xci. One can easily check that the maximum energy 7max ~ C jx„ ~ a^'^'^(.x-cr/.x'F)''"^^^^^ corresponds to (89). As 
demonstrated in Fig. 2., the analytical estimate (89) is in good agreement with the numerical integration of equation (93). 



5.5.5 Physical interpretation of particle acceleration 

At first glance, it is quite unexpected that a drag force can result in an acceleration of particles (so even the word 'drag' itself 
is not appropriate any more). As was demonstrated, additional acceleration can be realized for a highly magnetized outflow 
when the photon density decreases rapidly with radius (n > 2). This acceleration occurs as a result of the action of the drag 
force onto magnetic surfaces. Thus, acceleration is only obtained when the disturbance of the magnetic surfaces is taken into 
account self-consistently. 

To understand the nature of the additional particle acceleration in the supersonic region r > tf, it is necessary to return 
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to the system (90)-(91). As one can see, equation (91) actually plays the role of the GS equation, describing the force balance 

in the direction perpendicular to the magnetic surfaces. It contains no drag term, because the drag force acts along the 
magnetic surfaces. Indeed, as was shown before (Bogovalov 1998; Okamoto 1999; Beskin & Okamoto 2000), in the asymptotic 
region r ^ rp the transfield GS equation can be written down as !Fc = J'em, i.e., as a competition between the centrifugal 
volume force 

and the electromagnetic volume force 

= PcEo + -hB^ ~ ^. ^ ^ [{Bl - El) sin^ O] . (98) 
c 87rr sm Oti '- -' 

Here S « (c/A-k)EsB^ is the Poynting flux. Now using the expression for the curvature radius -Rc (Begelman & Li 1994; 
Beskin & Okamoto 2000), 

(99) 

Rc r sm or \ or J 
and the condition 

bI - El « ^bI (100) 

resulting from the relativistic Bernoulli equation (see e.g. Bogovalov 1998), for a magnetically-dominated flow S ^ nmc^'y we 
recover relation (91). 

On the other hand, the Bernoulli equation (90) describes the change of the total energy flux (45) due to the drag force. 

As was shown above, for Za > la in the vicinity of the fast magnctosonic surface the leading terms in the energy equation 
(90) are those containing ef and Za. This means that here the drag force again acts mainly on the electromagnetic (Poynting) 
flux S. The drag diminishes the ^-component of the electric field Ee, i.e. it disturbs the equipotential surfaces 6{r, 6) = const. 
But since in the one-fluid approximation the magnetic surfaces ef(r,6) = const must follow the equipotential surfaces, the 
decrease of the electric field Eg results in a change of the curvature of magnetic field lines as well. As the condition J^^ = T^m. 
can now be rewritten in the form 

(101, 

we see that the incresise of the curvature radius 7?c faster than r results in an increase of the particle energy. In our case such 
an acceleration is due to the rapid decrease (n > 2) of the isotropic photon density with the distance r. 

Thus, we see that for high enough photon density there is an additional acceleration of outgoing plasma outside the fast 
magnetosonic surface. On the other hand, as for Za > Zmax (64) the drag force significantly diminishes the total energy fiux 
E, one can conclude that this acceleration may only take place up to the energy 

1 ,o // s 4(™-2)/n(5-n) 

£sup ~ o^l^ (^) mec^ (102) 

for 2 < n < 3. On the other hand, fsup ~ (j^^^ttIoC^ for n > 3. This energy is always much lower than anieC? , corresponding 
to the total conversion of electromagnetic energy intoparticle energy. 



4 THE ELECTRON-PROTON OUTFLOW 



In this section we briefly consider the results of the analysis of the electron-proton outflow. As the procedure is quite similar 
to the electron-positron case, we present here the principal relations only. 

The main difference from the electron-positron case that occurs due to the large mass ratio is that the drag force acts on 
the electron component only, but the mass flow is determined entirely by protons. As a result, the energy conservation law 
has the form 



2e ^_ (Aj-_cose/2)7" -|-(mp/me) (A -cos 61/2)7+ 



tan^ 



Ia 



2\a sin 6 



(A+icose) (7-)'- 

Vme/ 



2Acr sin 

+ \2 



(A+icose) (7+) 



2e 
tan^"' 



7 



2a sin 6 2a sin Q 



Ia 



PX 

I u(x' 

J xa 



)^'^{x')Ax', 



(103) 
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max/^eC'^ a) positrons+electrons, n = 3 ^max/^eC'^ protons+electrons, n = 3 




Figure 3. Dependence of the maximum particle energy ^max on tlic compactness parameter Ip^ botli for electron-proton (ct) and electron- 
positron (6) outflows for n = 3. Here mp is a proton mass. The decrease of the particle energy for a very large compactness parameter 
/a is due to the decrease of the total energy flux . 



where now 7 = 7"*" . The conservation of angular momentum (47)-(48) looks like 



7+(l-xsin0^J) = 7+-Za 



u{x'){l - x sva.e$,+){-y+fAx - 4Acr— (5 - e/), 



7 (l-a;sin0$^) = 7i„ - 



/ u{x'){l - x' sine$-)(7-)"da;' + A\a{5 - ef), 

Jxn 



resulting in = l/a;sin^. Since for r rp we again have = ^,fi/{xsm6), we return to the universal dependence 
7 « a; sin 6. 

Equation (54) determining the energy on the fast magnetosonic surface now has the form 



3 o 'Tie 

7 — 2a 

rup 



2e/cos^-esin0|^-i^ 



i{x')j^{x')dx' + — 



2a;2 



2 , me .2/1 r. 

7 H asm 9 = Q. 

rUp 



Hence, for arbitrary Ia the Lorentz factor of particles can be presented in the form 

1/3 

7f 



(104) 
(105) 

(106) 

(107) 

(108) 



Now, for Ia < liv^ the position of the fast magnetosonic surface and the disturbance of magnetic surfaces on the fast 
magnetosonic surface are 



Xf 



2me 



2mc 



1/3 



1/3 • -1/3 
a ' sm ' 



-2/3 



-2/3 



Here for n < 3 we have 



;(P) 



™ N (5-n)/3 



\meJ 



.(n-2)/3 



(109) 
(110) 

(111) 



and for n > 3 we have — {mp/mc)^^^a^^^ . On the other hand, for Ia > Zcr we return to the relations xf = (ct/Za)^^^ 



(£/)f = lA/a for n > 3, but for n < 3 now 

, l/{5-n) 



Xf 



(112) 



(113) 



As we see, for the disturbance of magnetic surfaces and the position of the fast magnetosonic surface remain the same as 
for the electron-positron outflow. In particular, the terminating compactness parameter Zmax = cr is again determined by 
electrons. 

Thus, one can conclude that inside the fast magnetosonic surface r < rp the Lorentz factors of all particles are given 
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by the universal relation (106). Hence, our results demonstrate that the drag force does not affect the structure of the flow 

inside the fast magnetosonic surface for both the electron-positron and the electron-proton cases. But the position of the fast 
magnetosonic surface and the critical value of the photon density depend on the mass of the outgoing particles, and for the 
electron-proton case they are determined by protons. 

At large distances r > rp for Ia < Icr ^ the particle energy actually remains the same as on the fast magnetosonic surface 
{■y{r > rp) ?a 7p), while for Ia > & the Lorentz factor increases up to the majcimum value 

/ X 4(n-2)/(5n-n2) 

7max « 7F iyj^j (114) 

for n < 3. This takes place up to the distance r « rcr, where 

Clearly, this additional acceleration can be realized only for a > m^/me. Then, we have for n = 3 



^sup ' 



(T ' meC . (116) 

\meJ 



5 DISCUSSION AND ASTROPHYSICAL APPLICATIONS 

Wo have demonstrated how for a simple geometry it is possible to determine a small radiation-drag-forcc correction to the 
one-fluid ideal MHD outflow. The disturbance of magnetic surfaces was self-consistently taken into consideration. As a result, 
it is possible to characterize the general influence of the drag action on the magnetic field structure for an ideal magnetically- 
dominated quasi-monopole cold outflow and to determine under what circumstances radiation drag is important. 

As demonstrated above, the characteristics of the flow are determined by two main parameters, namely the compactness 
parameter I a (42) (which is proportional to the photon density) and the magnetization parameter a (34). If the photon density 
is low, so that the compactness parameter is small I a ^ ^cr(o"), the action of the drag force is negligible, while for a high 
photon density S> lci{(j), particles are additionally accelerated outside the fast magnetosonic surface. 

In particular, for Icr -C /a -C /max the increase of the drag force results in an increase of the outgoing plasma energy 
fmax ~ 7max»Tio.pC'^, but the disturbance of magnetic surfaces is small (e/ ^ 1). For Ia ~ 'max an increase of the photon 
density results in the increase of collimation up to values e/ ~ 1, but the particle energy remains near the saturation value 
fsup. Finally, for a very high photon density I a ^ /max an effective collimation of magnetic surfaces becomes possible, but 
in this case the drag force substantially diminishes the flux of electromagnetic energy inside the fast magnetosonic surface. 
As a result, for Ia 2> /max almost all the energy of the electromagnetic fleld is lost via the inverse Compton interaction of 
particles with external photons. For this reason, the very existence of a magnetically-dominated flow becomes impossible. The 
dependence of the maximum particle energy fmax ~ 7maxme,pC^ on the compactness parameter Ia is shown in Fig. 3. 

We now consider several astrophysical applications. 



5.1 Active Galactic Nuclei 

For AGNs (the central engine is assumed to be a rotating black hole with mass M ~ 10® Mq, R ~ 10^* cm, the total luminosity 
L ~ 10*^ erg s~^. Bo ~ lO'' G) the compactness parameter I a (42) can be evaluated as 

/a «30M9-i (^^) L45. (117) 
In the Michel magnetization parameter a (34) 

a « 10^*A-^M9B4 (^) , (118) 

the main uncertainty comes from the multiplication parameter A, i.e., in the particle number density n. Indeed, for an 
electron-positron outflow this value depends on the efficiency of pair creation in the magnetosphere of a black hole, which is 
still undetermined. In particular, this process depends on the density and energies of the photons in the immediate vicinity of 
the black hole. As a result, if the hard-photon density is not high, then the multiplication parameter is small (A ^ 10 — 100; 
Beskin, Istomin & Pariev 1992; Hirotani & Okamoto 1998). In this case for {Q.R/c) ~ 0.1-0.01 we have a ~ 10^ - 10^^ so 
that Zcr ~ 10^-10''. On the other hand, if the density of photons with energies E-^ > IMcV is high enough, direct particle 
creation 7 -|- 7 — > e"*" -f- e~ results in an increase of the particle density (Svensson 1984). This gives a ~ 10 — 10^, and hence 
Icr ~ 10 for an electron-positron outflow. 
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Prom a theoretical point of view, the most interesting result here is the possibility of an additional acceleration of particles 
outside the fast magnetosonic surface. Indeed, for a high enough photon density (I a ~ 10 - 100, i.e., for L ~ 10'"' - 
erg s~^) and a small magnetization parameter cr 10 — 100, the compactness parameter can exceed the critical value 
lev for an electron-positron outflow. In this case, according to Fig. 16, our analysis suggests that the kinetic luminosity of 

2/3 

the relativistic jet should be proportional to IJ^ oc Ltot, where Ltot in the total luminosity of the central engine. Kinetic 
luminosity is not easily determined from observations. However, observational evidence suggests that the radio luminosity of 
the jets is positively correlated with the luminosity of the central engine and the scatter of this correlation decreases towards 
larger luminosities (Baum, Zirbel & O'Dea 1995). 

For an electron-proton outflow the magnetization parameter (34) can be rewritten in the form (Camenzind 1990) 

. = f ^ « 3 X 10^ m)\lMl r \ " . (119) 

\ c ) cM \ c J ^ ^ VO.lMo/yry ^ 

Here M = 4nnmpR^c is the mass ejection rate. Hence, for a high ejection rate (M > 0.1 yr ~^) the magnetization 

parameter a < mp/me. In this case there is no acceleration of plasma. On the other hand, for low ejection rate M < 0.1 Mq/yt 
the magnetization parameter becomes too large for the drag force to be efficient. 

Thus, the drag force can substantially disturb the MHD parameters of a Poynting-dominated outflow only for a very 
high luminosity of the central engine (Ltot ^ 10*^ erg s ~^) and only for an electron-positron outflow. In all other cases the 
action of the drag force remains negligible. In particular, the additional acceleration of particles outside the fast magnetosonic 
surface is not efficient. 



5.2 Cosmological Gamma-Ray Bursts 

For cosmological gamma-ray bursts (the central engine is represented by the merger of very rapidly orbiting neutron stars or 
black holes with M ~ M©, R ~ 10'' cm, total luminosity L ~ 10®^ erg s"\ Bo ~ 10^^ G; see, e.g., Lee et al 2000 for details) 
the compactness parameter is extremely large: 

U ~ 10^^ (^) Lb2. (120) 

On the other hand, even for a superstrong magnetic field of Bo ~ 10^^ G (which is necessary to explain the total energy release) 
the magnetization parameter a is small (tr < 1 — 10), because within this model the magnetic field itself is secondary and its 
energy density cannot exceed the plasma energy density. Thus, one can conclude that for these characteristics of cosmological 
gamma-ray bursts the density of photons is very high so that Ia ^ Zmax and the drag force can make it difficult to form a 
Poynting-dominated outflow. A self-consistent analysis should include into consideration other physical processes such as high 
optical thickness resulting in the diminishing of the photon density, radiation and particle pressure, etc. Nevertheless, in our 
opinion, our conclusion may substantially restrict some recent models of cosmological gammas-ray bursts. 



5.3 Radio Pulsars 

For radio pulsars the central engine is a rotating neutron star with A4 ~ Af©, R ~ 10^ cm, total luminosity of the surface 
Lx - 10^^-10^'' erg s"\ and Bo ~ 10^^ G. In this case the magnetization parameter a ^ 10*-10^, corresponding to relativistic 
electron-positron plasma, is known with rather high accuracy (see, e.g., Bogovalov 1997). This gives Zcr ~ 10^-10^, and the 
compactness parameter 

Ia ~ (^) L35 (121) 

remains small (< 1) even for the most energetic (Lx ~ 10^^ erg s ~^ ) fast {flR/c ~ 10~^) pulsars like Crab and Vela. Thus, 

one can conclude that the drag force does not substantially disturb the magnetically-dominated outflow from radio pulsars. 

Thus, the drag force does not affect the wind characteristics (particle energy, magnetic field structure, etc.) of pulsars. 
However, interaction of outflowing relativistic particles with thermal photons can be important in other ways. In the wind 
region (r ^ _Rl) even a weak interaction with photons can result in a detectable flux of inverse Compton gamma-ray photons 
(Bogovalov & Aharonian 2000). On the other hand, near the surface of the star (r <ti Rh), inverse Compton photons are 
important in the pair creation process (Kardashev, Mitrofanov & Novikov 1984; Zhang & Harding 2000). 
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